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Abstract 

We have implemented a new numerical method to obtain the low-energy 
many-particle states of the Coulomb glass. First, this method creates an 
initial set of low-energy states by a hybrid of local search and simulated an- 
nealing approaches. Then, systematically investigating the surroundings of 
the states found, this set is completed. The transition rates between these 
states are calculated. The connectivity of the corresponding graph is analysed 
in dependence on temperature and duration of measurement. We study how 
the formation of clusters is reflected in the specific heat as non-ergodic effects. 
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Disordered systems of interacting localized particles have been extensively studied for 
over two decades [^0. A characteristic feature of these systems is a complex many valley 
structure of the energy landscape of the state space. In semiconductor physics, the Coulomb 
glass is a prominent example. It plays an important role as a semiclassical model for a 
disordered system of localized states with negligible quantum tunneling between them. 

We consider a half-filled impurity band in the strongly localized regime represented by 
the standard tight-binding Coulomb glass Hamiltonian 



where G {0, 1} denotes the occupation number of site i. The values of the random poten- 
tial ej are uniformly distributed between —W/2 and W/2. is the distance between sites i 
and j according to periodic boundary conditions. The sites are arranged at random with a 
minimum separation between them, which we choose to be 0.5 r, where r = (47rp/3)~^/^, and 
p is the concentration of sites. We take r as unit of distance, and the Coulomb interaction 
over a distance 1 as unit of energy. 

It is a comphcated task to obtain the low-energy many-particle states [0-0] • We have 
implemented a new numerical method which comprises sophisticated local search , ther- 
mal cycling ||^, and a renormalization approach to combinatorial optimization 0. First, 
we quench states chosen at random by means of a local search procedure, ensuring stabil- 
ity with respect to excitations on one up to four sites. Thus an initial set of metastable 
states is created. It is improved by cyclically heating (performing a fixed number of suc- 
cessful Metropolis steps), and quenching 0. In the course of this process, the temperature 
is decreased stepwise. Finally, we complete the set of low-energy states by systematically 
investigating the surroundings of the states found 

The set of low-energy states obtained is the basis for the study of different low- 
temperature properties. As a check, we verified that our results for the equilibrium specific 
heat agree with 0,0]. However, our main aim is to study the influence of the duration of 
measurement, r^, on the speciflc heat values, c. For that, we have to analyse the rates of 
the transitions between the states. 

The corresponding transition time (inverse of the transition rate in equilibrium) between 
two many-particle states / and J is a product of an energy factor and a spatial factor |Q, 



where k is the Boltzmann constant, taken as 1, and T the temperature. Ejj = max(£'7, Ej), 
where the ground state energy is assumed to be 0. The sum concerns only the sites, which 
change their occupation in the transition; it is the minimized sum of the related hopping 
lengths, a denotes the localization radius, and tq is a constant of the order of the inverse 
phonon frequency, tq ~ 10~^^ s. 

We consider the graph of the transitions between the many-particle states: Connec- 
tions indicate that r/j < r^. Thus the states are grouped into clusters. Assuming that 
thermalisation has happened inside the clusters, we measure the speciflc heat of a cluster a, 
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where the sum is performed over all states in a, and {...)a means thermal average in this 
cluster. 

We obtain the total specific heat, c(T, Tm), as weighted average of the c^. The weight Pa 
depends on the experimental situation simulated. We consider two situations: 

A) The 'sample' is at equilibrium. Thus Pa = Za/Z where Z denotes the partition 
function. Fig. |l| shows c(T = 0.012, r^) in comparison to the equilibrium value of c, where 
size effect, and reliability region are illustrated. Its main result is that also for durations of 
1 s to several hours the specific heat is significantly smaller than the equilibrium value. 

B) The 'sample' has been quenched from infinite T to the measuring T within a short 
time interval. To simulate this we quench first to T = 0, and heat then immediately to the 
measuring T: We start assigning the same probability to all states, and connecting them if 
the relaxation time (eq. without energy factor) is lower than the quenching time Tq. Note 
that this graph differs from the 'equilibrium graph'. 

Starting with the highest state of the cluster considered, we distribute its weight ac- 
cording to the transition probabilities to the states of lower energy. This process continues 
iteratively until only the local minima have a finite occupation probability. Finally, we assign 
to each 'equilibrium cluster' the sum of the probabilities of the included 'non-equilibrium 
local minima'. 

In Fig. we show the comparisons of both 'sample' preparation scenarios. Considering 
'samples' prepared as described in the previous paragraph, we obtain almost the same result 
as for 'samples' being in equilibrium (case A). Thus the question of the choice of the Pa is 
not an important task for the and Tq considered. 

We would like to acknowledge financial support from the SMWK, and from Acciones 
Integradas HA96-0065. 
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FIGURES 



FIG. 1. Dependence of the specific heat, c, on the duration of the measurement, Tm, for three 

sizes: O, □, and Q denote 64, 216, and 512 sites, respectively. The broken Hne represents the 
equihbrium value. The localization radius is a = 0.1 r, W = 2, T = 0.012, and tq 10~^^ s. 
Typical error bars are represented for one point of each curve; the ensemble averaging took into 
account 200 'samples'. 

FIG. 2. Comparison of two experiments with different initial conditions: Q = equili- 
brated 'samples', □ = quenched 'samples'. Here a = 0.1 r, 1^ = 2, 216 sites, T = 0.018, and 
Tg = 10^^ To = 0.1 s. Typical error bars are represented for one point of each curve; the ensemble 
averaging took into account 200 'samples'. 
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